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Abstract 

A microscopic model of financial markets is considered, consisting of many 
interacting agents (spins) with global coupling and discrete-time thermal bath 
dynamics, similar to random Ising systems. The interactions between agents 
change randomly in time. In the thermodynamic limit the obtained time series 
of price returns show chaotic bursts resulting from the emergence of attractor 
bubbling or on-off intermittency, resembling the empirical financial time series 
with volatility clustering. For a proper choice of the model parameters the prob- 
ability distributions of returns exhibit power-law tails with scaling exponents 
close to the empirical ones. 
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In recent years, both the analysis and modeling of financial time series have attracted a 
growing interest in statistical physics [1-11], motivated by a quest for a possibly universal 
dynamics underlying different markets. Indeed, various financial time series share similar 
properties, e.g., relative price changes (returns) show intermittent occurrence of large bursts 
(volatility clustering), resulting in the power-law scaling of tails in their probability distri- 
butions, as well as of the autocorrelation function of their absolute values ||. Large-scale 
microscopic simulations and analytic studies [4-11] revealed that such scaling can result 
from phenomena as diverse as, e.g., percolation |4| or stochastic processes with multiplica- 
tive noise ||[],|T^]. It was also shown P,|7|,|9|JTI1|| that in some models of market dynamics 
volatility clustering resembles on-off intermittency, an extreme kind of intermittency char- 
acterized by a sequence of large chaotic bursts separated by almost quiescent laminar phases 
|13| , or attractor bubbling, the predecessor of on-off intermittency in the presence of additive 
noise 
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In this Letter we present simulations of financial time series based on a microscopic model 
of many interacting agents with discrete-time dynamics. The agents are treated as spins 
exposed to fluctuations (a "heat bath") and coupled by randomly time-dependent Ising-like 
interactions. Related spin models, but with time-independent interactions, have been pro- 



posed for economic [TT| and social systems, where they are theoretically and empirically 
supported by the social impact theory of opinion formation |]L5|| . In contrast, we emphasize 
that the social interactions and the communication network among individuals are time 
dependent. They are also more and more determined by different kinds of long-range com- 
munication instead of spatial neighborhood which makes the topological structure of the 
interaction network unimportant. We show that, if the average strength of these interac- 
tions (which reflect the average reaction of agents to price changes) varies randomly in time, 
the system displays volatility clustering, and probability distributions of returns resemble 
those of empirical financial time series. A peculiar feature of our model is that the dynamics 
on a microscopic level is stochastic, but interactions between large numbers of agents cause 
the macroscopic dynamics to become typical of dynamical systems with attractor bubbling. 
In particular, in the mean-field approximation our model reduces to a generic model for 
attractor bubbling. Another important feature is the presence of various sources of stochas- 
ticity in the microscopic dynamics. First, the agents make decisions (on investments) under 
the influence of the external environment and other agents. The variation of the strength 
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of this influence is simulated by explicit random updates of the interaction strengths with 
the environment and among agents at every time step. Besides, the final agent decision 
is determined according to a probabilistic rule that mimics the the uncertainty in human 
decision-making. 

Explicitly, we consider i — 1,2, ... N agents (spins) with orientations o~i(t) = ±1, corre- 
sponding to the decision to sell (—1) or to buy (+1) a share of a traded stock or commodity 
at discrete time steps t. The orientation of agent i at time t + 1 depends on the local field 

1 N 

h{t) = -Y.M t )°A t ) + Ht), (i) 

where (t) are time-dependent interaction strengths among agents (reflecting their social or 
communication network), and hi(t) is an external field reflecting the effect of environment 
(e.g., access to external information which can differ between agents). The interaction 
strengths and external fields change randomly in time, 

Aij(t) = A£{t) + a Vij {t) , (2) 

Kit) = hQ(t), 

where, for simplicity, we assume that £(£), and Q(t) are random variables with no 

correlation in space and time, uniformly distributed in the interval (—1,1). According to Eq. 
(§), A is a measure of the randomly varying average interaction strength between agents, 
a is a measure of the random deviations of the individual interaction strengths Aij(t) from 
the average, and h is a measure of the random influence of the environment. The strengths 
Aij (t) can assume both positive and negative values, corresponding to a tendency to imitate 
or avoid the orientation of the interaction partners in the market. Below we show that the 
noises £(£) and Q(t) are necessary to get results comparable to empirical findings on price 
returns, while the effect of i]ij(t) is averaged to zero for large N. The orientations of all 
agents are updated synchronously according to the probabilistic rule 

{1 with probability p 
(3) 
— 1 with probability 1— p, 

describing an uncertainty in decision-making, where p — 1/ {1 + exp [— 2J»(t)]} (analogously 
to the well-known thermal bath dynamics). 
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In the following we will be interested in modeling the time series and distributions of 
the returns G^tit) = In S(t) — In S(t — At), i.e., forward changes in the logarithms of prices 
S(t) over a time scale At. Relative price changes are proportional to the difference between 
demand and supply, i.e., between the number of buying and selling decisions. Therefore, 
after introducing the average orientation x(t) = iV _1 J2iLi o~i(t) of agents we get dS/dt oc xS, 
which after discretizing time yields 
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G A t(t) oc Yl x(t-r) . (4) 

T=0 

In particular, Gi(t) oc x(t), which enables us to interpret the terms in Eqs. (0) and (p|) as 
follows: Since the average interaction strength A£(t) is common to all connections, it is a 
measure of the average reaction of the agents to the price changes, influencing their decisions 
via the term A£(t)x(t) in the mean field Ii(t). The terms arjij(t) describe the fluctuating 
interaction network while hQ(t) describe the fluctuating environment. The local field in 
Eq. (HD describes how agent % anticipates price changes by averaging the information on the 
opinions of his interaction partners and the external information accessible to him. 

Before discussing the results of numerical simulations for large N, we study the mean- 
field approximation of our model. For this purpose let us consider the case A ^ 0, a = 
and Q(t) = ((t) for all i. Then the dynamics of x(t) reduces to a one-dimensional map, 

x(t + l) = tanh [Af (t)x(i) + h((t)] (5) 
« A£(t)x(t) + Kit) , 

where the approximate equality holds for \A£(t)x(t)+ h((t) \ 1. The map (§) is a generic 
model for on-off intermittency if h = ||13|| , and for attractor bubbling if h > [I4|; let us 
also note that the approximate linearized form of Eq. @ belongs to a more general class 
of stochastic processes with multiplicative noise [^,|1^,[T(J] . In the case h = and with £(t) 



distributed uniformly in the interval (—1,1), the fixed point x = of Eq. (Bp loses stability 
for A > e = 2.718 . . ., as a result of the so-called blowout bifurcation, while for h > the 
stability is lost already for A > 1 |1~3| , ^4| . Above the instability threshold the trajectory 



departs, on average, from zero, but, due to the confining nonlinearity of the map (|5D, after 
the departure is randomly reinjected towards x = 0. Thus the time series of x(t) have 
intermittent character and consist of laminar phases, during which x(t) oc 0(h), and chaotic 
bursts, during which x(t) oc 0(1). Although the time series for on-off intermitency and 



attractor bubbling look similar, there is a qualitative difference between these phenomena: 
In the case of on-off intermittency the bursts are caused solely by the multiplicative noise £(t) 
whose amplitude A exceeds the blowout bifurcation threshold, while in the case of attractor 
bubbling the multiplicative noise only amplifies to macroscopic sizes small deviations from 
the fixed point caused by the additive noise C(t). In Fig. 1(a), an example of time series 
x(t) is shown in the attractor bubbling regime. It bears certain resemblance to empirical 
financial time series of price returns: There are frequent large bursts, corresponding to 
volatility clustering, separated by long-lasting laminar phases during which the dynamics of 
the map is governed mainly by the additive noise h((t). 

Another property of the map is that the probability distribution p(x) of x (pro- 
portional to the probability distribution of G\) shows pronounced tails obeying power-law 
scaling. It was shown analytically that for h = the distribution of x following the lin- 
earized multiplicative stochastic process in discrete time ([5]) obeys the power law scaling 
p(x) (x [fL6|J ; this scaling remains valid also for x ^> h > at least for Gaussian noises 

£, i] and continuous time |16|,[I7| . Following Ref. |i| the exponent a can be obtained as a 
non-zero solution of equation Jf^exp (ay) w(y)dy = A a (1 + a) 1 = 1, where y = ki\A£\ 
is a random variable with density w(y) = A~ 1 e y in the range y < In A. This means that 
the cumulative distribution P(x) exhibits power-law tails with P(x) oc x~ a . The power 
law scaling of p(x) and the analytic values of a (independent of h) are confirmed by our 
numerical simulations (Fig. 1(b)). In order to obtain a value a > 2 (i.e., within the range 
observed in empirical stock market time series H]), the parameter A in Eq. (||) should be 
in the attractor bubbling regime 1 < \A\ < e, so additive noise with h ^ is necessary. 
Moreover, the cumulative distributions of normalized returns G& t show tails obeying power- 
law scaling with robust exponent a for time scales At varying over more than one order of 
magnitude (Fig. 1(c)), a fact also typical of financial time series 0. 

It is also known that, close to the blowout bifurcation, the power spectrum of \x(t)\ 
from Eq. (]5|) shows power law scaling P(to) oc cu^ 1 ^ 2 for small uj, suggesting a similar power 
law scaling of the autocorrelation function of the absolute values of returns. However, in 
the interesting attractor bubbling regime we observed only slow exponential decay of the 
autocorrelation of \x(t)\. 

Now, we turn to numerical multi-agent simulations of the system. Examples of time 
series of the mean orientation x(t) of agents are shown in Fig. 2(a). It can be seen that, as 
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the number of agents increases, laminar phases and chaotic bursts become more distinct, in 
contrast with many other multi-agent models of stock market ||. In our model, the laminar 
phases with x ~ correspond to the disordered phase with no preferred orientation, and the 
bursts with a; oc 0(1) correspond to the ordered phase with a majority of agents sharing one 
orientation ("herding behavior"). For large N, and even for a > A, the time series of x(t) 
resulting from the simulation of Eqs. (1) to (3) and of the mean-field Eq. ([5]) look similar. 
Also the scaling exponents a are close to those obtained from the mean field dynamics (|Sp 
with the same value of A. Note that, in a system with a large but finite number of agents, 
there are always fluctuations of x(t) around zero in the disordered phase due to the heat- 
bath dynamics of individual agents. These fluctuations play the role of additive noise in Eq. 
(H), enabling behavior typical of attractor bubbling even if h = in Eq. (0). 

The main outcome of our many-body simulation is that, in a system of many agents, 
each of which is subject to "thermal bath" fluctuations, attractor bubbling can appear as 
a collective phenomenon, as a result of time- dependent interactions among agents. The 
requirement for its appearance is the average interaction strength between agents to vary 
randomly in time around the value at which the phase transition between the disordered 
and ordered state occurs. In the thermodynamic limit, the macroscopic system dynamics 
reduces then to the mean-field dynamics governed by the average interaction strength. The 
details of the interactions among pairs of agents, and the details of the variation of the mean 



interaction strength fl3| , seem unimportant for the qualitative result. 

To summarize, in this Letter we generalized a model of social impact to the case of 
randomly-in-time varying interactions among agents, and we have applied it to simulations 
of financial time series. Certain empirical properties of these time series, like volatility clus- 
tering and a power law distribution of returns, were shown to result from attractor bubbling, 
which appears in our model in the limit of many interacting agents. Thus we perceive the 
market as a dynamical system with an extreme kind of intermittency, in contrast with other 
models with time-dependent interactions between agents ||, which are related rather to 
self-organized criticality. Due to the lack of topological structure of the interaction network, 
the mean field approximation of our model quantitatively reflects the multi-agent dynam- 
ics. The necessary conditions for the applicability of the model are (i) random fluctuations 
in time of the average interaction strength between agents (corresponding to their average 
reaction to price changes), (ii) an uncertainty of decision making analogous to thermal heat 
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bath dynamics, and (iii), in the thermodynamic limit, small additive noise simulating the 
effect of the external environment. 

The authors would like to thank the Quandt Foundation of the ALTANA AG for financial 
support of the project "Nonlinear Dynamics in Models of Complex Systems". This work 
was supported in part by the Polsih Committee for Scientific Research under grant No. 5 
P03B 007 21. 

Fig. 1 (a) Time series of returns x(t) from the map (0) with A = 2, h = 10~ 2 ; (b) Probability 

distributions p(x) cx p{G\) of returns from the map (0) and fits to the power law scaling p(x) cx 

x-°- 1 of their tails for h = 1CT 3 and A = 1.6 (o, a = 2.67, a an = 2.89); h = 1(T 2 and A = 1.6 (□, 

a = 2.62, a an = 2.89), A = 2.0 (A, a = 0.98, a an = 1.0), A = 2.75 (O, a = -0.10, a an = -0.032); 

a an denotes analytic values of a. (c) Cumulative distributions of normalized returns G&t (where 

the angular brackets denote time average) from the map (|) with A = 1.6, h = 10~ 2 , At = 1 (o), 

At = 10 (□), At = 100 (A), At = 1000 (O), and comparision with the power law P (G A t) oc G~£ 

with exponent a = 2.76. All distributions were obtained from positive returns. 

Fig. 2 (a) Time series of returns x(t) from the model with many interacting agents with A = 2.0, 
a = 4.0 and, from above, N = 25 and iV = 4000; (b) Probability distributions p(x) of returns and 
fits to the power law scaling p(x) oc x~ a ~ l of their tails from the model with iV = 1000, h = 0, 
a = 2A, and A = 1.6 (□, a = 2.59), A = 2.0 (A, a = 0.98), A = 2.75 (O, a = -0.13). All 
distributions were obtained from positive returns. 
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